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The Millennium Galaxy Catalogue: The M^— L sp heroid 
derived supermassive black hole mass function 
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ABSTRACT 

Supermassive black hole mass estimates are derived for 1743 galaxies from the Millen- 
nium Galaxy Catalogue using the recently revised empirical relation between super- 
massive black hole mass and the luminosity of the host spheroid. The MGC spheroid 
luminosities are based on R 1 /™ -bulge plus exponential-disc decompositions. The ma- 
jority of black hole masses reside between 10 6 Mq and an upper limit of 2 x 10 9 Mq. 
Using previously determined space density weights, we derive the SMBH mass func- 
tion which we fit with a Schechter-like function. Integrating the black hole mass func- 
tion over 10 6 < Mbh/M Q < 10 10 gives a supermassive black hole mass density of 
(3.8 ± 0.6) x 10 5 h% M Q Mpc~ 3 for early-type galaxies and (0.96 ± 0.2) x 10 5 tf M Q 
Mpc~ 3 for late-type galaxies. The errors are estimated from Monte Carlo simula- 
tions which include the uncertainties in the M^-L relation, the luminosity of the 
host spheroid and the intrinsic scatter of the M^-L relation. Assuming supermassive 
black holes form via baryonic accretion we find that (0.008 ± 0.002) /i 3 per cent of the 
Universe's baryons are currently locked up in supermassive black holes. This result is 
consistent with our previous estimate based on the M\fa—n (Sersic index) relation. 
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1 INTRODUCTION 

Observations of nearby galaxies have shown that super- 
massive black holes (SMBHs) exist within the centres of 
both ellip tical galaxies and the classic a l bulges of disc 
galax ies l|Kormendv fc Richstonel 1 19951 : iMagorrian et al.1 
1 19981 ). At present there are about 50 credible measure- 
ments a n d 26 partial measurements as summarised in 
I Graham! (|2008bh . For lower mass systems (< 10 s Mq), 
galaxy core s are dominated by a luminous compact mas- 
sive object (iGraham fc Guzm an 2003; Ferrarese et al.|[2006l ; 



iBalcells et al.l l2007f ) . iGallo et al.1 (|2008h found AGN activ- 
ity in early-type galaxies with M ga i < 10 10 Mq , revealing 
that SMBHs c a n exi st in low mass galaxies. Furthermore, 
ISatvapal etall (|2008h found SMBHs can form and evolve 
in bulgeless galaxies. These argue t hat SMBHs have built 
up from less massive 'seeds' (e.g., IWoo et all l2008r i. The 
SMBH se eds are believed to fo rm from either PoP III stars 
collapse (IVolonteri et al.l 120031) o r directly from dense gas 
collapse (jKoushiappas et al T 12004I : iBegel man et al.ll2006l ) 111 
which large quantities of gas are driven to the galaxy's core 
giving rise to an accreting black hole or active galactic nu- 
cleus (AGN). Eventually, the energy outflow from the AGN 



may become sufficient to halt further gas infall and truncate 
further star formation in the surrounding stellar spheroid, 
either thro ugh bright /je t mod e or r adio/diffuse mode feed- 
back (e.g. lCroton et al.1 120061 and ISomervillel l200Sf ). Al- 
though some pape rs al s o adv ocate AGN activity as a trigger 
for star formation Isiikl (120051 ) and lPipino et al.1 (|2009l ). This 
interplay between the SMBH and its surroundings is be- 
lieved to give rise to correlations between the SMBH mass 
and a variety of measurable properties of the host galaxy's 
spheroid component^ After initial formation, SMBH growth 
is believed to progress through periodic gas accretion onto 



1 Some properties known to correl ate with SMBH mass 
are: the spheroid l u minos ity L jKormendv fc Richstonel 

ersion 



1 19951 : iMarconi fc Huntl l2003t): the mean vel ocity disp' 



a of the spheroid ijFerrarese fc Merrittl |2000| : iGebhardt et al.l 
120001: iGrahaml |2008al; iHul | 20081) : the spheroid ma ss M aD h 
I KormendT'fcRi^hstoiie] Il995l : IMarconi fc Hunt! l2003t 
lHaring fc Rbd |2004|) :~ the galaxy light concentration C 
jGraham" eTaUl20*0l|): the Sersic index n of th e majo r-axis surface 
bright ness profile llGraham fc Drived |2007|) . See iNovak et al.1 
1120061 ) for a comparative review. 
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Figure 1. A comparison of SMBH masses derived from the M^—n relation versus the M^-L relation (left panels: z < 0.06, m iddle 
panels: z < 0.1, right panels: z < 0.18). The red circles are the MGC data. The blue circles denote spheroids from lCaon et all \ l99 3T ). The 
spheroids less luminous than Mg = —18 mag or with Sersic index < 2 (indicated with dashed vertical and horizontal lines, respectively) 
are shown as grey squares as these data are less reliable. The sloped dashed line indicates the 1-1 relation. Top panels: (left to right) 
30 - 210 - 564 elliptical galaxies (B/T = 1), middle panel: 39 - 312 - 867 early- type disc galaxies (1 > B/T > 0.4), bottom panel: 34 - 
109 - 312 late-type disc galaxies (0.4 > B/T > 0.01). The "right" error bars in each panel show the average uncertainties on the masses 
derived from of the intrinsic scatters of the M^-L,n relations and the "left" errors shows those derived from the errors on the spheroid 
luminosities and Sersic index (see Section 3.1). The correlation coefficient (r) has been estimated for each subsample. 



the SMBH during major me r ger events (e.g. , Granato et alj 
|2004 |Pi Matteo et al.ll2005l ; Hopkins et alll2008l ). This re- 
sults in the coalescence of the SMBHs and spheroids in 
such a manner as to preserve these init ial correlations to 
the present time (|Robertson et al1l2006l ). Recent evidence 
does however suggest some evolution in the SMBH scaling 
relations and so weakening t he criterion that these relations 
must be perfectly preserv ed (|Shields et al.|[200?3 ; IWoo et all 
l200Sl ; IShankar et all [20091 ). 

AGN activity peaks in the range 1 < z < 2 (|Ueda et al.l 
120031 ; iRichards etalll2006l ; iHopkins et Id] 120071 ) which sets 
the key formation epoch for massive bulge/spheroid forma- 
tion and the key merger phase of massive dark matter haloes. 
If the picture above is correct then the zero-redshift SMBH 
mass function and the c osmic AGN en ergy history should be 
closely coupled. In fact. lSoltanl l|l982h argued that the AGN 
cosmic history and local SMBH mass function can be com- 
bined to determine the energy efficiency of the accretion pro- 
cess. Recent attempts in this direction have been made by a 
number of groups (e.g.. iTamura et all 120061; I Hopkins et al 



l2006l , 120081 ; IShankar et all 120091 ; ICao fc Lil l2008l ; lYu fc~ 



l200ct ) leading to significant conclusions. However, the cred- 
ibility of these studies rests heavily on the robustness to 
which the local SMBH mass function has been measured 
and to which the cosmic AGN history is known. 

Over recent years there have been a number of esti- 
mates of the SMBH mass function at zero or low redshift 
and a c ompendium of some recent results are presented in 
Fig. 5 oflShankar et all (|2009T ). Except for our previous work 
(|Graham et al.ll2007h none of these studies actually empiri- 
cally measure the SMBH mass function but instead estimate 
it by analytically combining the measurements of the galaxy 
luminosity or velocity function with on e of the SMBH scal- 
ing relations (as concisely outlin ed bvlHaxing fc Rixll2004l . 
following the first application by ISalucci et al. 19991 ). This 
inherently incorporates some assumptions of, for example, 
the bulge-to-disc ratio and how it varies with luminosity or 
how galaxy luminosity relates to galaxy velocity dispersion. 

In this paper we depart from this 'analytic' approach 
and pursue a direct empirical strategy wh ich uses t he re- 
cently revised M^-L relation provided by iGrahaml (|2007l , 
his equations 6 & 19) to obtain individual SMBH masses for 
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each galaxy in t he Millennium Galaxy Catalogue (MGC; 
iLiske et al.ll2003l ) which has a reliable spheroid luminosity. 
This complement s our previous empir ical method using the 
Mbh-n relation (|Graham et al.l l2007h and we discuss how 
the two results compare. 

In Section[2]we describe the data and sample selection. 
In Section [3] we review the latest M^-L relation. In Sec- 
tion [4] we use this relation to calculate the SMBH mass for 
each galaxy and find the SMBH mass function of local galax- 
ies, the mass density of SMBHs, the cosmological SMBH 
density Obh and the baryon fraction of the Universe that is 
contained in central SMBHs. Finally, in Section [5] we com- 
pare our results to those previously derived from the MGC 
using the M^ - n rela tion and those recently summarised by 
IShankar et ail (|2009l ). 

A cosmological model with Qa = 0.7, fi m = 0.3 and 
Ho = 70 ft 70 km s -1 Mpc -1 is adopted throughout this 
paper. 



2 THE MILLENNIUM GALAXY CATALOGUE 



The Millennium Galaxy Catalogue fMGC: lLiske et al.ll2003l ) 
is a medium-deep survey (fiu m = 26 mag arcsec -2 ) cover- 
ing a wide region of sky (37.5 deg 2 ) in the B-band (4407 A). 

The survey extends 75 deg along the equator (from 9 h 58 m to 
14 h 47 m) The data frames 

were obtained using the 4-CCD 
mosaic Wide Field Camera on the 2.5-m Isaac Newton Tele- 
scope. Each CCD has a pixel scale of 0.333 arcsec pixel -1 . 
The data were taken with a median seei ng FWHM = 1.3". 

Using the SExtractor package |Bertin fe Arnoutsl 

1 19961 ) a catalogue was derived containing over one million 
objects in the range 16 ^ -Bmgc < 24 mag. Initially, two 
catalogues were defined: MGC-BRIGHT which includes all 
galaxies with B < 20 mag and MGC-FAINT containing 
the rest. For fu rther details on th e MGC imaging data and 
its analysis see ILiske et"al] (j2003l ). MGC-BRIGHT is com- 
prised of 10095 resolved g alaxies and has 96 per cent com- 
plete redshift information (|Driver et al.l [20051 ). This sample 
was decomposed into bulges and discs with GIM2D using 
an Ft 1 / 71 Sersic profile for b ulges and an exponential profile 
for discs (| Allen et al.ll2006T ). The robustness of the decom- 
positi on process has be en verified using duplicate observa- 
tions (lAllen et al"1l2006f ) and by using extensive simulations 
(| Cameron et al.l 20091 ). All data and data products used in 
this paper are freely available^ 



2.1 Sample selection 

For this study we use the online catalogue mgc-gim2d and 
extract the following parameters: spheroid absolute magni- 
tude, bulge-to-total (B/T) flux ratio, half-light radius of the 
bulge, redshift and weight (see Section 4). We convert the 
absolute magnitudes from Ho = 100 km s -1 Mpc -1 to 70 km 
s _1 Mpc - . We additionally restrict our sample to the red- 
shift range 0.013 to 0.18. As this sample may contain galax- 
ies with nuclear components (e.g. star clusters) that bias the 
R 1/n model, we remove all 'spheroids' with half-light radii 



http:/ /www. eso.org/~jliskc/mgc/ 



and B /T less than 0.333 a rcsec (1 pixel) and 0.01, respec- 
tively (|Graham et al.ll2007l ). We note that the MGC region 
is overdense by 9.08 per cent (Hill et al. 2009, in prep.) and 
we normalise our weights accordingly. 

The sample we con struct here is identical to sample 3 
of lGraham et al] (|2007l ) which incorporates a colour cut of 
(u — r) coro > 2.0. This colour cut isolates the red spheroid 
population and is intended to remove the predominantly 
blue pseudo-bulge systems. Pseudo-bulges are believed to 
form by a process distinct to that by which classical bulges 
form, and it is not yet clear whether they contain a SMBH 
in their centre and, if so, whether they ad here to th e pre- 
viously cited spheroid-SMBH relations (see Irlull2008l). The 
colour cut adopted here follows the findings in Driver et al.l 
(|2007h of a clear colo ur bimodality within sp heroidal systems 
(as also identified bv lDrorv fc Fisherll2007l ). After we apply 
the above colour cut we divide the remaining galaxies into 
two subsamples. The first contains 1431 'early-type' galax- 
ies (B/T > 0.4) and the second contains 312 'late- type' 
galaxies (0.01 < B/T < 0M. 

Finally, we note that I Allen et al.l (|2006l ) consider bulges 
with M b > — 17 mag to be less reliable while iGraham et al.l 
|2007l ) remove all galaxies fainter than Mb = — 18 mag from 
their sample. In this paper we indicate the 'unsafe' area with 
a dashed vertical line at log(Mbh/M0) = 7.67 in all relevant 
figures which corresponds to Mb = — 18 mag. Data fainter 
than this limit should be treated with caution. 



3 THE Mbh-L RELATION 

The first empirical correlation between the mass of the 
SMBH a nd the bulge luminosity was ide ntified in the review 
paper bv lKormendv fc Richstonel (|l995l ). Following this ini- 
tial study a lot of similar relations have been derived but 
with differing slopes for the Mbh~L relation. As a result, for 
the same bulge luminosity we have a range of possible SMBH 
masses. This large degree of scatter that has appeared in the 
luminosity relation is exacerbated by the difficulties in esti- 
mat ing the b u lge lu minosity in mid-type galaxies. 

iGrahaml (|2007h reviewed the relations of four studies 
and updated their samples with new data (including re- 
vised estimates of the distances, Hubble type classifications 
and SMBH masses), added new galaxies, and removed some 
which still have significant uncertainty in their listed values. 
In this paper we use t he best fit from a 22 galaxy B-band 
sample (|Grahamll2007l . his equation 19): 

log(Af b h/M Q ) = -0.40(±0.05)(M s +19.5)+8.27(±0.08), (1) 

with a total scatter of 0.34 dex in log Mbh and an intrin- 
sic scatt e r of .30 dex. These 22 galaxies were derived by 
Graham! ||2007|) fr om a parent sample of 27 galaxies from 
Marconi fc Hunt! (|2003l ) Group I galaxies. It is important 
to note that this relation does not contain any corrections 
for internal dust attenuation. However, recently we demon- 
strated that dust attenuation can amount to as much as a 
few m agnitudes in B depending on inclination (|Driver et al.l 
l200ct ) . This poses somewhat of a dilemma as we have a choice 
of whether to use the final relationship presented by Graham 
and our original uncorrected MGC magnitudes and space 
densities (as used in our previous estimate of the SMBH 
mass function) or whether to use the dust-free Mhh~L re- 
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Table 1. SMBH mass function data for the full, early- and late- type samples (DC = dust-corrected), as presented in Fig. [2] The 
uncertainties given are the la values derived from the Monte Carlo simulations. 

log 10 M bh [10- 4 h 3 7(j Mpc" 3 dex" 1 ] 

[Mq] All galaxies Early-type Late-type All galaxies (DC) Early-type (DC) Late-type (DC) 
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lation for elliptical galaxies (|Grahamll2007l . his equation 6) 
combined with our dust corrected MGC bulge magnitudes 
and the correspondingly revised space densities. We choose 
to use and show the results of both alternatives in order to 
indicate the possible uncertainty introduced by dust atten- 
uation. 



3.1 Masses from Mbh— L versus masses from Mbh— n 

iGraham et al. I (|2007l ) derived the SMBH mass of each galaxy 
in the MGC sample using the photometric quantity Sersic 
index0 (n). In this paper we derive the SMBH masses using 
the spheroid luminosity. In Fig. [JJ we directly compare the 
SMBH masses derived for each galaxy using the two inde- 
pendent methods (Mbh~£ and Mhh-n) for various subsam- 
ples defined in B/T and redshift as indicated. A correlation 
between the two BH mass estimates is not seen and there- 
fore requires explaining. For each subsample we measure the 
linear correlation coefficient r and only the low-z elliptical 
systems show evidence for a correlation. For the bulges of 
disc galaxies, no correlation is detected. 

However, we find that the lack of correlation can be well 
explained when one considers the errors from both the mea- 
surements (on L and n), the intrinsic scatters of the Mbh~ 
L,n relations, and the limited mass range probed here (see 
error indicators on Fig[T}. In more detail, equation (TTJ above 
implies that t he typical bulge l uminosity error of ±0.1 mag 
(see fig. 15 of lAllen et al.ll200fj ) corresponds to a mass error 
of ±0.04 dex. The internal scatter of the M^h~L, however, 
is ±0.30 dex. The Sersic index value has a typical measure- 
ment error of ±20% for the elliptical sample and ±35% for 
the early & late- type sam ple. While the intern al scatter of 
the M hh -n is ±0.18 dex (|Graham et al.l 120071 '). Hence the 
scatter in Fig. [T] is simply representative of the combination 
of two measurement errors and two intrinsic scatters (with 
the measurement error dominating for the Mbh — n estimates 



3 For more infor mation on the Sersic index including a compre- 
hensive review see IGraham fc Driverl l l2005h . 



and the intrinsic scatter dominating for the Mbh — L esti- 
mate) coupled with the relatively narrow range of parameter 
space probed. 

To illustrate this conclusion we have used the same 
Mbh~L,n relations as above to calculate the SMBH masses 
for a sample of nearby, lo w-inclination syst ems within the 
Virgo and Fornax clusters (|Caon et al.l|l993l ), and for which 
the measurement errors are, by comparison, negligible. Over 
the full range of parameter space explored, in this study, 
these data are known to give relatively tight correlations. 
However when overplotted (blue points on Fig. [TJ for the 
narrow range of parameter space we probe here, they give 
similarly low correlation coefficients as the MGC data. This 
strongly suggests that the lack of correlation is mainly due 
to the limited range of parameter space covered by the MGC 
data. The lack of correlation is consistent with the adopted 
errors (and which are propagated throughout the analysis). 
The caveat is that the analysis remains susceptible to any 
unknown systematic errors (e.g., bar-bulge contamination) 
however the overlap between the MGC data and local data 
suggests these systematic errors cannot be a dominant fac- 
tor. Finally we note that the dominant error in using the 
Mbh — L relation comes from the inherent intrinsic scatter 
in the M b h — L relation rather than the MGC measurement 
error. Conversely the reverse is true in our previous study 
based on the Mbh — n relation. 



4 THE MGC SMBH MASS FUNCTION VIA 

Mbh - L 

We have derived individual black hole masses for each 
spheroid using equation JTJ. These estimates have three 
sources of error. The first is the uncertainty on the param- 
eters defining equation (UJ). This is a systematic error. The 
second source of error is the uncertainty in our spheroid 
magnitude measurements. Finally, as noted in the previous 
section, the M b h~i relation has significant internal scatter. 
We consider the latter two as random errors. (Note that the 
alignment of the calibration sample with the MGC sample 
on Fig. [TJ implies that there is no obvious systematic error 
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Table 2. Results of a three-parameter Schechter function fit to the SMBH mass function, both for dust corrected (DC) and uncorrected 
samples. The uncertainties on the parameters were derived from Monte Carlo simulations (see text for details). 



Sample 


No. of 


lo£ 


4>- 




log(M,/M ) 




a 






spheroids 


[^70 M P C 


-3 


dcx^ 1 ] 














All galaxies (B/T > 0.01) (DC) 


1743 


-3.15 


± 


0.05 


8.71 


± 


0.08 


-1.20 


± 


0.10 


All galaxies (B/T > 0.01) 


1743 


-2.92 


± 


0.05 


8.71 


± 


0.06 


-1.07 


± 


0.05 


Early-type (B/T > 0.4) (DC) 


1431 


-3.10 


± 


0.05 


8.64 


± 


0.03 


-0.90 


± 


0.06 


Early-type (B/T > 0.4) 


1431 


-2.90 


± 


0.04 


8.71 


± 


0.03 


-1.16 


± 


0.10 


Late-type (0.01 < B/T < 0.4) (DC) 


312 


-3.70 


± 


0.02 


8.46 


± 


0.03 


-1.11 


± 


0.11 


Late-type (0.01 < B/T < 0.4) 


312 


-3.66 


± 


0.02 


8.49 


± 


0.08 


-1.28 


± 


0.01 
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Figure 2. The data points in the top panel show the SMBH mass 
function of the full MGC sample. The error bars were derived 
from Monte Carlo simulations (see text for details). The solid 
line and grey shaded region represent the best fit Schechter func- 
tion and its uncertainty. The dashed line shows the same for dust 
corrected m agnitudes and weigh ts, while the dotted line shows 
the result of lGraham et all ||2007V ) using the M^—n relation (cor- 
rected for the 9.08 per cent overdensity of the MGC). The middle 
and bottom panels show the same for our early- (B/T > 0.4) and 
late-type (0.01 < B/T < 0.4) samples. The dashed vertical line 
at log(Mt,h/A4"0) = 7.67 (corresponding to Mb = —18 mag) indi- 
cates our reliability limit and we show all data beyond this limit 
with big open circles. The fits use only the reliable data above 
this limit. All data points are list ed in Table [T] The a dditional 
curve in the middle panels is from lHopkins et al 1 (|2007l 'l (H07). 
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p bh (h? o lO»M Mpc- 3 ) 

Figure 3. A histogram of the 10001 values of mass density for 
early-type (B/T > 0.4) galaxies from the Monte Carlo simula- 
tions. The distribution is reasonably reproduced by a Gaussian 
(solid line). 

in our flux measurements.) To model these errors appro- 
priately we run a series of 10001 Monte Carlo simulations. 
For each simulation we randomly modify the coefficients of 
equation fl]) assuming Gaussian error distributions. We then 
perturb each spheroid flux by a random amount as is appro- 
priate given the spheroid's individual magnitude error, and 
calculate a new set of SMBH masses from the perturbed 
fluxes using the perturbed version of equation (JTJ) - Finally, 
we randomly modify each of the SMBH masses according 
to the internal scatter of the Mbh-L relation, again assum- 
ing a Gaussian distribution. Given these 10001 Monte Carlo 
datasets we estimate the value and la error of any given 
quantity of interest by calculating that quantity for each of 
the datasets and determining the median and 68 percentile 
range of the resulting distribution, respectively. This pr oce- 
dure is similar to that followed bv lGraham et all (|2007| ). 

We now wish to volume-correct our sample and to com- 
pute the SMBH mass function. T o this end we emplo y the 
space-density weights derived by iDriver et al.l (|2007l ): the 
weight of a spheroid of luminosity L is simply given by the 
value of the luminosity function of the appropriate spheroid 
type, 4>(L), divided by the observed number of spheroids in 
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the interval [L, L + dL]. The value of the SMBH mass func- 
tion at mass Mbh is then calculated as the sum of the weights 
of all spheroids within the interval [Mbh, Mhh + dMbh]: 



0(M bh ) = 



E 

[M,M+dM] 



W(L), where W(L) = <j>(L)/N(L). (2) 



In Fig. [5] we show the SMBH mass function and its 
errors as determined from the Monte Carlo process described 
above and using equation[2]for our full, early- and late- types 
subsamples (top to bottom). All data are listed in Table [T] 
In each panel of Fig. [2] the solid line represents the best fit 
of a three-parameter (M» , 0* , a) Schechter-like function over 
the mass range 10 7 75 < M bh /M < 10 10 : 



<KM bh ) = 0* 



exp 



Mbh 
M* 



(3) 



The best fitting parameters are given in Table [2] 

Finally, we integrate our Schechter-like estimation of 
the SMBH mass function over the mass range probed by 
our data to calculate the SMBH mass density, given by the 
expression: 

/•log(M bh /M Q) =10 

Pbh = / ^(Mbh)MbhdlogAfbh 

Jlog(M hh /M @) =6 



^M.e 1 
ln(10) 



7 + 2, 



10 1U M« 



AL 



7 U + 2, 



10 6 A/ G 

A'L 



(1) 



As an example of the Monte Carlo process, Fig. [3] shows 
the distribution of the 10001 SMBH density values. This 
is reasonably described by a Gaussian function. The final 
SMBH density is defined as the median of this distribu- 
tion with the la error given by 68 percentile range. The 
resulting mass densities of SMBHs with masses between 
10 6 < Mbh/M© < 10 10 for the full, early- and late-type 
samples are given in Table [3] The last column of Table [3] 
also lists the corresponding cosmological density parameters 

fibh = pbh/ pcrit- 

Integrating the mass function over all masses we derive 
slightly larger values for the densities: 5.06 (8.5), 3.8 (6.6) 
and 0.96 (0.92) xlO 5 ft? M Q Mpc" 3 for the full, early- and 
late-type samples, respectively, where the numbers in paren- 
theses refer to the dust-corrected valu es. These densities ar e 
consistent with the values reported bv lGraham et all (|2007h . 

As discu s sed in the introduction and also in 
IShankar et all lj2004) if one assumes that SMBHs form via 
the accretion of baryons alone and using the cosmological 
SMBH mass density fibh.totai = ( 3 -7 ± 0.7) x 10~ 6 h 7Q we 
can estimate the SMBH baryon fraction. Given that 4.5 
per cent of the critic al density is in the form of baryons 
(|Tegmark et alj|2006h we find that (0.008 ± 0.002) h% per 
cent of the Universe's baryons are currently locked up in 
SM BHs at the centres of g alaxies. For comparison, the value 
that lGraham et alj (|20071 ) found was (0.007 ±0.003) h% per 
cent. 



5 DISCUSSION AND COMPARISONS TO 
EARLIER WORKS 

Fig. U shows a compendium of recent estimates of the zero 
redshift SMBH mass function in the left panel, while the 
right hand panel shows the differential contribution to the 



mass density. The shaded (light grey) region representing 
the MGC results show the spread in the means from our 
three estimates as t abulated in Table [5] and in Table 1 of 
IGraham et~ail l|2007l ). Within this compendium of data only 
the MGC results are based on actual measurements of bulge 
properties. The non-MGC curves/shading all rely on cou- 
pling an empirical SMBH mass-relation with a galaxy lu- 
minosity or velocity function under some assumption as to 
the fraction of ellipticals and how the mean bulge-to-total 
ratio varies with luminosity. Typically a constant fraction 
for both values is adopted independent of luminosity which 
is contrary to what has been seen in the MGC and in other 
morphological studies. In reality the fraction of ellipticals 
and the mean B/T ratio should increase with stellar mass. 
Incorporating this would have the effect of tilting the mass 
functions derived in this 'analytical' way, and this is the 
most likely explanation for the discrepancy between some 
of the results shown in Fig. U and those from the MGC 
(see Appendix A for deta ils). This echoes the r ecent study of 
nearby active systems bv lGreene fcHol \200H ) who also find 
a bounded BH mass function based on an empirical study of 
the broad-line regions of 9000 QSOs, albeit shifted to lower 
black hole masses (see Fig. |4j left). 

However, it is also fair to note that while the MGC re- 
sults give a consistent SMBH mass density and generally 
agree at the high mass end they show a broader discrep- 
ancy at lower SMBH masses. Most likely this can be traced 
back to the bulge-disc decompositions which are no doubt 
imperfect, particularly for low-B/T, low-luminosity and/or 
low-n systems (i.e., those systems likely to contain lower 
mass SMBHs). Repea tability tests (Ullen et all 120061 ) and 
extensive simulations (I Cameron et alJl2009h based on MGC 
data both show that while disc parameters are recovered 
extremely reliably the bulge parameters are susceptible to 
gross error. However, no obvious bias via the simulations 
has yet been seen, but rather a larger than desirable general 
scatter and in this manner the correct answer with appro- 
priate errors should indeed emerge through the Monte Carlo 
simulation process. Note that the modelling of the SMBH 
mass function includes extensive Monte Carlo simulations 
of the structural errors along with a Malmquist Bias correc- 
tion. 

Within the MGC results the dominant source of ran- 
dom error comes from the uncertainty in the Mbh~n and 
Mbh~L relations and the associated intrinsic scatter. This is 
limited mainly by the sample size of the relevant calibration 
datasets. In terms of systematic error this most likely arises 
from an, as yet unknown, bias in the bulge-disc decomposi- 
tions. While simulations show that systematic biases are low 
we do see a discrepancy between the two approaches which 
at this stage provides the best and only indication of the 
scale of any systematic. Simulations suggest that our bulge 
luminosity measurements are more robust than our Sersic in- 
dex measurements (reflected by the lower final errors in the 
SMBH mass function) and so we therefore consider the new 
M GC measureme nts based on the updated Mbh~L relation 
of lGrahaml (|2007l ). to be superior to our previous estimate 
and to the other 'analytical' estimates (given the inherent, 
sometimes hidden assumptions in the analytic approach). 
We therefore advocate the data presented in this paper as 
the current best estimate of the zero redshift SMBH mass 
function. 
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Table 3. SMBH mass densities integrated over the mass range 10 6 - 10 10 for the full, ear ly- and late-type samples as derived from the 
M b ^—L (dust corrected [DC] and uncorrected) and M b ^—n relations llGraham et al,|[2007l) . The differences in the numbers of spheroids 
are due to different absolute magnitude cuts employed, a is the error due to the uncertainties on the parameters defining equation JT) 
and j3 is from the uncertainty on the bulge magnitudes. 7 is derived from the MGC global cosmic variance of 6 per cent for the effective 
30.8 dog 2 region with 0.013 < z < 0.18. Finally, <5 is the error due to the intrinsic scatters of the M b ^-L (0.30 dex) and M b ^-n relations 
(0.18 dex). All densities have been corrected for the overdensity of the MGC region (9.08 per cent). 



Method No. of 




Pbh ± a 


±/3±-y±5 






^bh 


spheroids 




[10 5 h 3 70 


M Q Mpc- 3 






[10- 6 /l 70 ] 


Full sample (B/T > 0.01) 














M hh -L 1743 


4.9 


± 0.7 ± 


0.5 ± 0.3 


± 


0.1 


3.7 ±0.7 


M hh -L DC 1743 


8.0 


± 1.3 ± 


0.4 ± 0.5 


± 


0.3 


6.2 ± 1.1 


Mbh-n 1543 


1.0 


± 1.5 ± 


0.06 ± 0.2 


± 


0.04 


2.9 ± 1.1 


Early-type (B/T > 0.4): 














M bh -L 1431 


3.8 


± 0.6 ± 


0.4 ± 0.2 


± 


0.1 


2.5 ± 0.3 


M hh -L DC 1431 


6.5 


± 1.2 ± 


0.3 ± 0.4 


± 


0.3 


4.8 ±0.6 


M bh -n 1352 


3.1 


± 1.04 ± 


0.05 ± 0.2 


± 


0.03 


2.3 ±0.8 


Late-type (0.01 < B/T < 0.4) 












M bh -L 312 


0.96 


± 0.2 ± 


0.05 ± 0.06 


± 


0.07 


0.7 ±0.2 


M bh -L DC 312 


0.92 


± 0.1 ± 


0.04 ± 0.05 


± 


0.05 


0.7 ±0.1 


M bh -n 191 


0.86 


± 0.49 ± 


0.03 ± 0.05 


± 


0.02 


0.6 ±0.4 



SMBH mass functions derived from the M - a rela- 
tion also show an increase in density towards low black 
hole masses. However velocity dispersions for the lower 
mass systems are generally estimated rather than mea- 
sured with a number of hidden im plicit assumptions re- 
quir ed (e.g.lAller fc Richstonell2002l) The most recent stud- 
ies (|Lauer et al.l 120071 ; iTundo et all I2007T) u se the veloc- 
ity dispersion function from lSheth et all 02OO3I . however this 
function estimates the black hole masses from the systemic 
circular velocity which could lead to an overestimation of 
SMBH masses and a steeper faint-end. 



6 CONCLUSIONS 

We have used a sample of 1743 galaxies extracted from the 
Millennium Galaxy Catalogue t o estimate the S MBH masses 
from the M bb -L relation (using lGrahamll200 71 . his equation 
19) to construct the SMBH mass function at redshift zero. 
Our results agree well with our previous study, however, in 
comparison to other published data the MGC data consis- 
tently show a steep drop-off or decline in the space-density 
of SMBHs towards lower masses. Although this lies below 
what we consider to be our reliable limits (< 1O 7 M0) it 
most likely arises from our more direct empirical approach. 
Essentially, the SMBH mass function declines because t he 
spheroid luminosity function declines (|Driver et al.l [20071 ) . 

Within our reliability limits the contribution of the 
SMBH mass function to the total cosmic SMBH density 
budget is sharply peaked. Integrating over the best fit mass 
function we obtain a total SMBH mass density of pbh = 
(4.90 ± 0.7) x 10 5 h 70 M Q Mpc" 3 and a cosmological SMBH 
density of Q bh = (3.7 ± 0.7) x 10" 6 h 70 . This implies that 
(0.008 ± 0.002) h 70 per cent of the Universe's baryons are 
contained in SMBHs, in excell ent agreement with the re- 
sults from iGraham et al.l (|2007l ). 

In this, and our previous, paper we have measured the 
SMBH mass fun ction via three met hods. Once using the 
Mbh-n relation (|Graham et al.l 1200 7| ) and twice using the 
Mbh-L relation (this paper; once ignoring dust attenuation 



and once correcting for dust attenuation). From a compari- 
son of the three SMBH mass functions we find that the dis- 
crepancy in any one bin of the SMBH mass function (prob- 
ably the most realistic quantitative insight into the hidden 
systematic errors) lies at the 30 per cent level. This consti- 
tutes the current limit which can be achieved from the MGC 
dataset. A significant improvement can made in three ways. 
Firstly, by increasing the sample size, secondly by obtaining 
more robust measurements, and thirdly by improving the 
accuracy in which the M b h-L,n relations are known. One ob- 
vious direction in the latter case is to re-derive the M^-Ljn 
relations in the NIR to overcome t he severe effects of dust 
attenuation (see lDriver et al.l|2007l ). At the present time no 
suitable catalogue exists. Over the coming years we are em- 
barking on a programme to recalibrate the M bb -L,n rela- 
tions at near-IR wavelengths (where we expect the intrinsic 
scatter to be significantly reduced), as well as construct a 
~ 100 x larger NIR-selected galaxy sample (~240k galaxies 
via the ongoing Galaxy And Mass Assembly survey; Driver 
2008). With these two improvements we should be able to 
significantly reduce the errors inherent in this work as well 
as probe a significantly broader black hole mass range. 
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Figure 4. Left panel: a com pendiu m of zero redshift SMBH mass functions co rrected for the hidde n Hubble de pendency. The curves 
arc from lShankar et al] l|2004h fS04): lAller fc Richstonel l|2002l ) (AR02a, AR02b~l: [Greene fc Hoi l|2007l l (GH07) and lshankar et al] l|2009h 
(S08,z=0 curve). The light grey shaded region represent the MGC resu l ts fro m our three estimates. The data points show the SMBH 
mass function for all the MGC subsamples (red squares: iGraham et al blue triangles: this study, dust corrected; green circles: 

this study, no dust corrected). Right panel: equivalent plot showing the contribution by SMBH mass to the cosmic SMBH mass density. 
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APPENDIX A: PREDICTING THE SMBH 
MASS FUNCTION DIRECTLY FROM THE 
GALAXY LUMINOSITY FUNCTION 

In this paper we have derived the nearby SMBH mass func- 
tion through direct estimates of SMBH masses for 1743 
individual galaxies selected from the Millennium Galaxy 
Catalogue (MGC; see Sample 3 definition described in 
iGraham etail l2007h . Our resulting SMBH mass function 
determined using the Mhh~L relation is consistent with our 
earlier estimate based on the Mhh-n relation and those de- 
rived analytically as shown in Fig. U However, below our 
limit of bulge reliability (i.e., Mb > — 18 mag) we note a 
sharp decline in our mass function not mirrored in the an- 
alytical estimates. In Section [5] we claim this is due to in- 
herent assumptions in the analytical process related to the 
accuracy to which type fractions are known and the adoption 
of a constant bulge-to-total luminosity ratio. To understand 
why this is so we firstly describe the analytical method, we 
then explore some of the implicit assumptions, and finally 
derive a set of analytical SMBH mass functions to illustrate 
the basis of these claims. We note that to fully explore this 
issue it is necessary to revert to the full MGC sample for 
which some fraction of our bulge measurements will be un- 
reliable (see lCameron et al.ll2009l for a detailed discussion of 
these issues). We therefore make a strong caveat that the 
values and results reported in this appendix are to quali- 
tatively illustrate some of the subtleties in estimating the 
SMBH mass function analytically only and are not suffi- 
ciently robust to produce definitive SMBH mass function 
estimates at this time. 



Al Estimating SMBH mass functions analytically 

Apart from our empirical studies, all estimates of the SMBH 
mass function combine a global galaxy luminosity function 
with estimates of population ratios and bulge-to-total lumi- 
nosity ratios. This situation arises because the Millennium 
Galaxy Catalogue is the only study, to date, to actually 
separate the bulge and disc c omponents for a sufficiently 
large sample (|Allen et al.ll2006l ). In the analytical approach 
using, for example, the M^-L method, a galaxy luminos- 
ity function is adopted (or a set of Hubble type specific lu- 
minosity functions) and the space-density of bulges is then 
estimated by adopting a constant bulge-to-total flux ratio 
(or a set of bulge-to-disc flux ratios for each Hubble type). 
These estimated bulge luminosity functions can then be con- 
volved with the Mbh~i expression to produce an analytic 
SMBH mass function. There are three main issues with this 
technique: (i) Hubble population fractions are not well con- 
strained and depend on luminosity (while type dependent 
luminosity functions will accommodate for this, estimates 
based on the global luminosity function typically do not), 
(ii) Bulge-to-total flux ratios are believed to vary with lumi- 
nosity (typically all systems fainter than Mb = — 17 mag are 
single component only) and this has not yet been adequately 
quantified, (iii) A final more subtle issue is the distinction 
between the type of bulge, simple Hubble classification does 
not distinguish between red bulges (classic) and blue bulge 
(possible pseudo- bulges or contaminant bars). Our assertion 
in this paper has been to only predict SMBH masses for red- 
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Figure Al. (main panel) The distribution of B/T versus abso- 
lute magnitude for various galaxy sub-populations (as indicated). 
Overlaid are the mean trends and errors on the mean for the red- 
bulge(solid squares), blue-bulge(solid triangles), and combined 
spiral population(solid circles), (upper) The relative fraction of 
various galaxy sub-populations colour coded as in the main panel. 



bulge systems as it is unclear at this time whether blue-bulge 
systems also harbor SMBHs. 



A2 MGC B/T dependency with luminosity 

Returning to the full publicly available MGC catalogue we 
can directly determine the dependency of bulge type and 
bulge-to-total ratio as a function of luminosity. First we re- 
strict the sample (to minimise the contamination by spu- 
rious bulges) to the redshift range 0.013 < z < 0.18 leav- 
ing 7044 galaxies from a parent sample of 10095. We now 
construct the following sub-populations (following the argu- 
ments laid down in earlier MGC papers): 

(1) Red bulge-only (elliptical) systems (B/T — 1, (u — r) c > 
2), 579 

(2) Blue bulge-only (elliptical) systems (B/T = 1. (it — r) c < 
2), 235 

(3) Red bulge systems (0 < B/T < 1, (it - r) c > 2), 2487 

(4) Blue bulge systems (0 < B/T < 1, (u - r) c < 2), 964 

(5) Bulge- less or unresolved bulge systems (B/T = 0), 2779 
We plot these populations on the main panel of Fig. IA1I 
Solid points show the mean B/T for the red and blue bulge 
systems as well as that for the combined spiral sample (red 
bulge, blue bulge and bulge-less systems) versus absolute 
magnitude. We see that the bulge-to-total fraction of the 
individual red and blue bulge populations are distinct but 
relatively constant. In contrast the overall mean spiral B/T 
shows a linear decrease (reasonably well described by the 
relation B/T = -1.09 - 0.067M B ). 

In the top panel of Fig. lAll we show how the fraction of 
each population varies. At bright luminosities the sample is 
dominated by red-bulge systems which steadily declines, in 



turn we see a corresponding rise in blue-bulge and bulge-less 
systems. The trends seen are quite distinct showing minimal 
scatter in the measurement of the means, this suggests that 
while errors in some of our bulge measurements are no doubt 
present they are unlikely to be dominating. 

A3 Derivation of the SMBH mass function 
through scaling relations 

We can now adopt a global galaxy lumin osity function (e.g. , 
that derived in B for the MGC itself in lDriver et al.1 120061 ) 
and combine this with the scaling relations seen in lAll under 
various assumptions as follows: 

Method 1 — Adopt a constant B/T and a constant spiral 
fraction (standard practice). 

Method 2 — Adopt a varying B/T and a constant spiral 
fraction. 

Method 3 — Adopt a constant B/T and the red bulge frac- 
tion (similar to our Sample 3). 

Method 4 — Adopt a constant B/T and the red and blue 
bulge fraction. 

Fig. IA2I (upper panel) shows the global MGC luminosity 
function (red/solid) and the implied bulge luminosity func- 
tions based on these four methods outlined above. Fig. IA2I 
(lower panel) shows the implied SMBH mass function. We 
can see that the faint end of the SMBH mass function de- 
pends critically on these relations and our assumptions. In 
particular the key question arises as to whether blue bulges 
system (potentially pseudo-bulge systems) harbor SMBHs 
or not and if so whether they follow the same trends as clas- 
sical (red) bulges. In conclusion, and having explored both 
methods, we continue to advocate the empirical approach 
laid out in the main body of this paper, where every bulge 
has been measured, and the errors have been derived from 
Monte Carlo simulations, over the analytical approach with 
its hidden assumptions. 
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Figure A2. (upper panel) The global B-band galaxy luminosity 
function from Driver et al. (2005) and a variety of bugle luminos- 
ity functions one might derive (as described in the text), (lower 
panel) The resulting analytical SMBH mass functions depending 
which model LF is adopted. 



